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Recent lattice QCD calculations have reported evidence for the existence of a bound state with 
strangeness —2 and baryon number 2 at quark masses somewhat higher than the physical values. 
By developing a description of the dependence of this binding energy on the up, down and strange 
quark masses that allows a controlled chiral extrapolation, we explore the hypothesis that this state 
is to be identified with the //-dibaryon. Taking as input the recent results of the HAL and NPLQCD 
Collaborations, we show that the H-dibaryon is likely to be unbound by 13± 14 MeV at the physical 
point. 
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Our understanding of quantum chromodynamics has 
been challenged for decades by the apparent absence of 
multi-quark states. The outstanding candidate for such 
a state has been the i?-dibaryon, ever since it was sug- 
gested that it should be very deeply bound with respect 
to the A — A threshold jT| . Extensive experimental efforts 
to find this new particle PHS] have led to the conclusion 
that it does not appear to be bound. However, the issue 
has been given new life in the past few months by reports 
from the HAL and NPLQCD Collaborations, whose lat- 
tice simulations find that the i7-particlc is indeed bound 
at quark masses somewhat above the physical range [Hl[7]. 

Early spectroscopic studies of the baryon and di- 
baryon spectrum within the bag model [S] showed that 
the pion cloud contribution was phenomenologically very 
important in determining whether or not the H was in- 
deed bound. A preliminary analysis of some of the lattice 
data using a fit linear in suggests that it is bound 
at the physical quark mass, while an extrapolation lin- 
ear in m7r, while also consistent with binding, does allow 
that the H may be slightly unbound [9]. We consider it 
important that the extrapolation to the physical quark 
masses should respect the constraints of chiral symmetry, 
such as ensuring the correct leading non-analytic behav- 
ior. Given the rather large range of quark masses over 
which current lattice simulations have been made, the 
technique which offers the best opportunity for a quan- 
titative fit, while preserving the correct non-analytic be- 
havior, is finite range regularization (FRR) [TDHH]. 

In this Letter we explore the possibility that the H is 
a compact, multi-quark state that may be bound with 
respect to the A — A threshold. We apply the FRR tech- 
nique to describe the quark-mass dependence of both 
the octet-baryon masses and the binding energy of the 
i/-dibaryon. With few lattice results available for the 
dibaryon binding energy, it is essential to utilise the fit to 
the hyperon masses to determine the dependence on the 
non-singlet combination of quark masses (i.e., mi — rris), 
which plays a critical role in describing the variation of 
the (non-singlet) A mass. Extrapolated to the physical 
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TABLE I. Values for the terms linear in the non-strange quark 
mass, mi — > m^/2, and the strange quark mass, — >■ {m\ — 
mJ/2), expressed in terms of the leading quark-mass insertion 
parameters q,/3 and a. 

quark masses, we conclude that the if-dibaryon is most 
likely unbound, with its mass being 13 ± 14 MeV above 
the A - A threshold. 

Following the technique described in Ref. [Uj, we 
fit the data for octet masses recently published by the 
PACS-CS Collaboration [TJ, using an expansion about 
the SU(3) limit for the light and strange quark masses: 

Mb = Af(o) + W4'' + <5Mf -I- . . . (1) 

Here the leading term, M^°' , denotes the degenerate mass 
of the baryon octet in the SU(3) chiral limit, and 

(5mW ^ - C^m,, (2) 

with the coefficients given in Table |Tj is the correction 
linear in the quark masses. 

At next order, after the linear mass insertions, one 
finds quantum corrections associated with the chiral 
loops involving the pseudo-Goldstone bosons, (f) = n, 
K, rj. While our formal assessment of the order of a 
given diagram treats the intermediate octet and decu- 
plet baryons as degenerate, in order to more accurately 
represent the branch structure near ~ (5, we retain 
the octet-decuplet mass difference (S) in the numerical 
evaluations. These loops take the form: 

Wff - E [xs0/fl("^0,O, A) 

+XT<i,lR{m^,S,A)] , (3) 
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where the coefficients xb4>, XTcft, taken from Ref. [TS], are 
given in Table [TTl 



The meson loops involve the integrals: 

Ir = - I dk = V? {k)-bo-h2m'i 

(4) 

where the subtraction constants, &o,2, are defined so that 
the parameters M(°\ C^^j and C^^l are renormalized (ex- 
plicit expressions may be found in Ref. [10 , or can be 
readily evaluated numerically by Taylor expanding the 
integrand in m^). The loop contribution parameters are 
taken to be £> + F .g^ = 1-27, = fZ?, C = -2D, 
f = 0.0871 GeV and S = 0.292 GeV. Within the frame- 
work of FRR, we introduce a mass scale, A, through a 
dipole regulator, u{k) — ( ^2^fe2 A is related to the 
scale (typically ~ A/3 for a dipole) beyond which a for- 
mal expansion in powers of the Goldstone boson masses 
breaks down [20[ I21j. However, rather than growing un- 
controllably as some power of m^, in this regime the 
Goldstone loops are actually suppressed, decreasing as 
powers of A/m^ [TBHT5] . In practice, we choose this 
mass parameter by fitting the lattice data itself. Exten- 
sive studies of the FRR technique have established that 
the extrapolation is independent of the functional form 
chosen for the regulator ID] - essentially because of the 
rapid decrease of the loop contributions which we just 
explained. 

In order to describe the mass of the H, treated as a 
compact, multi-quark state (rather than a loosely bound 
molecular state), we note that because it is an SU(3) sin- 
glet, at the equivalent order of the quark mass expansion 
its mass can be expressed as: 

Mh = m(,°' -aH{^+ml) + SM'i'^'^ , (5) 

with 

^^4^/^) = - Ch(2D2 _^ (^2) [/^(^^^ A) 

+ pB{mK, SH, A) + |/fl(m,„ 5H, A)] . (6) 

As indicated in Eq. because it is a flavour-singlet, 
at leading order in the quark masses depends only 
on the sum of the quark masses. We therefore set: 

Bh = 2Ma - Mh 

= (2Af (°) - a4°)) - (4aA - aH){'f + mj,) 

-2am% - 2pml + 2<5A/f - SM^"^^ 

2 

= So - '7b(^ + - ^aml - 2l3ml 

+25Mf'^^ - JA^/^^ , (7) 



where Bq and gb are parameters determined by the fit 
to the lattice data for Bh [S [7]. Of course, a and (i 
are determined by the fit to the baryon octet described 
above. For simplicity, we keep the regulator mass for 
the octet and the H the same, while varying the chiral 
coefhcient for the iJ, Ch, to fit the lattice data. 

In order to fully maintain the correlations between the 
errors associated with all of the fitting parameters, we 
carried out a simultaneous analysis by minimizing for 
the fit to both the masses of the nucleon octet and to 
the difference in mass of the H and two A hyperons. Of 
course, the parameters Af a, /3, a, as well as the regu- 
lator mass A, were primarily determined by the fit to the 
PACS-CS data, which is shown in Fig. [l] with the cor- 
responding parameters given in Table [Till We note that, 
as explained in Ref. [T^ , the octet data was corrected for 
small, model independent finite volume effects before fit- 
ting. In the figure we show the lattice data after applying 
both the finite volume correction and a correction (based 
upon our fit parameters) arising because the PACS-CS 
simulations used values of the strange quark mass that 
were somewhat larger than the empirical values. The 
per degree of freedom for the octet data was 0.49 (7.3 
divided by 20 — 5 = 15). This is lower than unity as, 
without access to the original data, we cannot incorpo- 
rate the effect of correlations between the lattice data. 
Nevertheless, the fit is clearly very satisfactory over the 
entire range of quark masses explored in the simulations 
and should provide an excellent basis for the study of the 
possible binding of the i?-dibaryon. Indeed, the masses 
of the N, A, S and S baryons at the physical point are 
(0.959±0.023 , 1.129±0.014 , 1.188±0.011 , 1.325±0.006) 
GeV, where all the errors include the correlated uncer- 
tainties of all the fit parameters, including the regulator 
mass, A. For comparison we note that the physical octet 
masses are (0.939 , 1.116 , 1.193 , 1.318) GeV. 



With respect to the i?-dibaryon we have retained only 
the data from the HAL Collaboration [B] which was gen- 
erated on the largest lattice volume, namely 3.87 fm. 
These data points correspond to large (degenerate) pseu- 
doscalar masses, 1.015, 0.837 and 0.673 GeV, for which 
the finite volume corrections are expected to be very 
small. Accordingly we used the reported values without 
applying any finite volume correction. In the case of the 
NPLQCD Collaboration [T, where the calculation was 
performed at TOtt = 389 MeV and tiik = 553 MeV, we 
include in the fit only the value for the binding of the H 
determined after their extrapolation to infinite volume. 

We chose the mass splitting between the ff-dibaryon 
and the other dibaryon states appearing in its chiral loop 
corrections to be the same as the octet-decuplet mass 
splitting used earlier, namely 6H = 5 = 0.292 GeV. This 
is compatible with the estimates of Aerts et al. [52] cal- 
culated within the MIT bag model, as well as with the 
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TABLE 11. Chiral SU(3) coefficients for the octet baryons to octet (B) and decuplet (T) baryons tfirough tfie pseudoscalar 
octet meson S. 
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TABLE IIL Values of the fit parameters for the octet and _ff-dibaryon data corresponding to the fits shown in Figs. 1 and 2. 




(GeV^) 



0.4 0.6 
m^- (GeV2) 



FIG. 1. Fit to the octet data of PACS-CS [13] using Eq. ([T]). 
Note that we have fit the data after applying finite volume 
corrections and we have also used our fit to correct the lattice 
data for the strange quark mass, which was somewhat larger 
than the physical value. 



FIG. 2. Binding energy of the B-dibaryon versus pion mass 
squared, resulting from our chiral fit, for several values of the 
kaon mass at which the simulations of Refs. [6] and [7] were 
carried out. 



experimental absence of other nearby states. The sen- 
sitivity of our fit to (577 is quite small, with an increase 
(decrease) of 8 by 100 MeV increasing (decreasing) the 
mass difference 2Mp^ ~ Mfj by only 4 MeV. This small 
shift is combined in quadrature with the error found from 
our chiral fit to yield the final, quoted error in the binding 
of the H. 

The best fit parameters describing the binding energy 



of the /7-dibaryon are also given in Table III and the ac- 
tual fit is shown in Fig. [2] As explained above, the data 
shown in Fig. [2] are from NPLQCD (lowest mass point) 
and HAL (three largest mass points). In each case the 
curve nearest the data point illustrates the extrapolation 
as a function of the light quark mass implied by our fit 
at the value of the strange quark mass corresponding to 
that lattice data point. The errors shown are the result 
of combining in quadrature the statistical and systematic 
errors quoted by the collaborations. The shaded error 



bands incorporate the effect of correlations between the 
fit parameters, including the uncertainty on the regula- 
tor mass. We note the remarkable result that the best 
fit value of the chiral coefficient for the iJ-dibaryon, Ch 
(which for convenience is normalized with respect to the 
chiral coefficient for tt loops on the A hyperon in Eq. ([6|), 
is within 20% of the theoretical value reported by Mul- 
ders and Thomas [8], who calculated it using SU(6) sym- 
metry. If instead we retain the Mulders-Thomas coeffi- 
cient, the i/-dibaryon is unbound by 30 ±9 MeV but the 
quality of the fit is significantly reduced, with a per 
degree of freedom of almost 2. 

It is clear from Figs. 1 and 2 that both the octet data 
and the data on the binding energy of the if-dibaryon 
are very well described. The binding of the H is reduced 
by a decrease in the masses of the u and d quarks and 
the s, with significant chiral curvature for ttItt below 0.4 
GeV. It is important to note that our analysis does not 
include the effect of the coupling of the _ff-dibaryon to 
the open A — A channel. Given where we find the multi- 
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quark state, this might be either attractive or repulsive, 
depending on the range of momenta which dominate the 
continuum channeL It is also true that it would be ex- 
tremely valuable to have new simulations at lower pion 
mass and large volume which would further constrain the 
extrapolation. Nevertheless, the conclusion of our study 
is quite clear. Even though the H is bound at larger 
quark masses, we see in Fig. [2] that chiral physics leads 
to a more rapid decrease of the mass of the A as ap- 
proaches its physical value than we find for the H and as 
a result one must conclude that at the physical values of 
the quark masses the iJ-dibaryon is most likely unbound. 
Our estimate, including the effect of correlations between 
all the fit parameters, is that the H is unbound by a mere 
13±14 MeV at the physical point. That this is so close to 
the A — A threshold will undoubtedly spur investigations 
into the consequences for doubly strange hypernuclei as 
well as the equation of state of dense matter. 
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